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Abstract. Theoretical spectral energy distributions computed for a grid of hydrogen-deficient and carbon-rich 
model atmospheres have been compared with the observed infrared (1-2.5 /xm) spectra of V4334 Sgr (Sakurai's 
Object) on 1997 April 21 and July 13. The comparison yields an effective temperature of T c g — 5500 ± 200 K 
for the April date and T c ff = 5250 ± 200 K for July. The observed spectra are well fitted by Asplund et al. 
(1999) abundances, except that the carbon abundance is higher by 0.3 dex. Hot dust produces significant excess 
continuum at the long wavelength ends of the 1997 spectra. 
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1. Introduction 

V4334 Sgr (Sakurai's Object), the "novalike object in 
Sagittarius" discovered by Y. Sakurai on February 20, 
1996 (Nakano et al. 1996) is a very rare example of ex- 
tremely fast evolution of a star during a very late final 
helium-burning event (Duerbeck & Benetti 1996). During 
the first few months after discovery, Sakurai's Object in- 
creased in visual brightness to V ~ 12 m . In 1997 it in- 
creased further to V ~ ll™ 1 . In March 1997 the first ev- 
idence of dust formation was seen ( Kimeswenger et al. 
1997, Kamath & Ashok 1999, Kerber et al. 2000). In 
early 1998 the optical brightness of Sakurai's Object de- 
creased ( dimming first reported by Liller et al. 1998), but 
then recovered. However, during the second half of 1998 
an avalanche-like growth of the dusty envelope occurred, 
causing a rapid decrease in optical brightness and the com- 
plete visual disappearance of the star in 1999. At present 
essentially only thermal emission by dust can be observed 
(Geballe et al. 2002). Our view of the born again star has 
been completely obscurred by the dust it has produced. 

Abundance analyses by Asplund et al. (1997, 1999) 
and Kipper & Klochkova (1997) have found peculiari- 
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ties similar to those of R CrB-like stars. Asplund et al. 
(1999) estimate that the logarithmic abundances of hy- 
drogen, helium and carbon in atmosphere of Sakurai's 
Object in October 1996 were -2.42, -0.02 and -1.62, re- 
spectively]], with hydrogen only the third most abundant 
element by number. All of the above studies are based 
on optical spectra obtained in 1996. At that time the 
spectrum of Sakurai's Object resembled that of an F- 
supergiant; molecular bands were absent or very weak. 
Cooling of the photosphere of Sakurai's Object resulted 
in its optical spectrum during 1997 and 1998 resembling 
those of C-giants with very strong bands of CN and C2 
(Pavlenko,Yakovina & Duerbeck 2000). Modeling of some 
of these optical spectra have allowed estimates of the 
changes in T e ff and Eb-v to be made during this period 
of rapid evolution of the optical spectrum (Pavlenko et al. 
2000; Pavlenko & Duerbeck 2001). 

Modeling of near infrared (1-2.5 /im) spectra of 
Sakurai's Object is of interest for several reasons. In addi- 
tion to providing comparisons with results obtained from 
the optical spectrum and tests of the reliability of molec- 
ular and atomic data, it allows accurate determination of 
the effective temperature and sensitive tests for emission 
by hot dust. Use of the 1-2.5 /im region for modeling is 
especially important after 1996, when the bulk of the pho- 
tospheric flux shifted from the optical into this waveband. 



In this work we will use an abundance scale ^2 Ni — 1. 
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In this paper we present and compare model 1-2.5 /im 
spectra with those of Sakurai's Object obtained during 
1997, on UT April 21 and July 13 at the United Kingdom 
Infrared Telescope (UKIRT). The observed spectra to- 
gether with observational details were presented by Eyres 
et al. (1998) and the July spectrum is also shown in 
Geballe et al. (2002). The resolutions of these spectra as 
presented here are 1.4 nm (0.0014 /im) at 1.02-1.35 /im 
and 2.8 nm (0.0028 /im) at 1.42-2.52 /im. Narrow spectral 
features in the 1.82-1.95 /im portions of these spectra are 
due to incomplete removal of strong telluric lines. 

2. Modeling Procedure 

Grids of plane-parallel model atmospheres in LTE, with 
no energy divergence were computed by the SAM12 pro- 
gram (Pavlenko 2002). This program is a modification of 
ATLAS 12 (Kurucz 1999). Opacities due to C I bound-free 
absorption, of importance in atmospheres of hydrogen- 
deficient, carbon-rich stars over a wide (0.1-8 /im) wave- 
length region (see Pavlenko 1999, 2002; Asplund et al. 
2000), were computed using the OPACITY PROJECT 
cross sections database (Seaton et al. 1992). The opacity of 
C _ also was taken into account Myerscough & McDowell 
(1996); see Pavlenko 1999 for more details). 

An opacity sampling approach (Sneden et al. 1976) 
was used to account for atomic and molecular line ab- 
sorption. The source of the atomic line information was 
the VALD database (Kupka et al. 1999). Lists of diatomic 
molecular lines of 12 CN, 13 CN, 12 C 2 , 13 C 2 , 12 C 13 C, 12 CO, 
13 CO, SiH, and MgH were taken from Kurucz (1993). We 
adopted Voigt profiles for every absorption line; damp- 
ing constants were computed following Unsold (1955). 
Microturbulent velocities of 3-6 km s _1 were adopted. 

Two grids of model atmospheres with different abun- 
dances were computed. One grid used the chemical com- 
position of Sakurai's Object determined by Asplund et 
al. (1999) for October 1996. The other is the same except 
that the abundance of carbon is increased by 0.6 dex. This 
"carbon-rich" case is of interest because the carbon abun- 
dance of Asplund et al. (1997) is larger by 0.6 dex than 
that which was obtained by them later from the analy- 
sis of high resolution spectra (see Asplund et al. 1997, 
1999, 2000. This "carbon problem" appears to arise more 
from the analysis of C I lines than from C II or C 2 lines 
(Asplund, private communication). For both grids the iso- 
topic ratio 12 C/ 13 C= 5 was adopted (Asplund et al. 1997). 

To determine molecular densities, a system of equa- 
tions of chemical equilibrium was solved for a mixture 
of ^70 atoms, ions and molecules, including the most 
abundant diatomic molecules containing carbon. We used 
the approach developed by Kurucz in ATLAS 12 (Kurucz 
1999) in which ratios of densities of atoms and n x , ... ,n z 
and molecules n xy ,_ z follow the equations: 

n x * ... * n z /n x ... z = exp(-D a /T ev + b- 
c*(T + d*(T-e*(T + /* T))) + 



3/2 * (to - fc - 1) * InT), (1) 

where Do and T ev are the dissociation potential and the 
temperature in eV (for neutral diatomic molecules m=l, 
k=0). The adopted molecular constants for the three most 
important molecules are given in Table 1. 

Synthetic spectra^ were computed using the WITA6 
program (Pavlenko 1999) for the same grid of opaci- 
ties, abundances, isotopic ratios and microturbulent ve- 
locites that were adopted for the model atmospheres. 
WITA6 computes spectra and spectral energy distribu- 
tions (SEDs) taking into account "line by line" absorption 
by atomic and molecular transitions. For Sakurai's Object 
spectra computed with wavelength step 0.05 nm were con- 
volved with gaussian profiles with full widths at half max- 
imum of 1.4 and 2.8 nm over the appropriate wavelength 
intervals. Computed and observed spectra were normal- 
ized at 1.7 /im for comparison; see Pavlenko, Yakovina & 
Duerbeck (2000) for more details. 

3. Results 

3.1. Principal spectral features 

In Fig. |l| the principal features formed by the molecular 
species C 2 , CO, and CN are displayed in separate spectra. 
Atomic features are also shown, as these are also present 
in Sakurai's Object (see Eyres et al. 1998; Geballe at al. 
2002). As in the optical spectrum (Pavlenko et al. 2000), 
absorption of only a few molecular species accounts for the 
main features in the IR spectrum. Only the most abun- 
dant isotopic species of each molecule is shown. Of the less 
abundant isotopic species, only bands of 13 CO have been 
detected in the infrared (Eyres et al. 1998). 

3.2. Dependences on T c g , log g, Vt and log N(H) 

The model spectra of Sakurai's Object display a strong de- 
pendence on T c ff (Fig. ||).|^]ln general, the dependence of 
the IR SED on T e g is determined mainly by the variations 
of the molecular densities with temperature. The band 
strengths of CN, CO and C 2 all increase as T c g decreases. 
Changes in the continuum fluxes are much smaller. Similar 
effects are seen in model optical spectra (Pavlenko & 
Yakovina 2000). However, there the molecular bands are 
numerous, whereas in the infrared only the few strongest 
vibration- rotation bands of CN, C 2 , and CO are promi- 
nent. 

As can be seen in Fig. ^, the dependence of the spec- 
trum on log g is generally rather weak. However, there are 
differences in the responses of different spectral regions. 
The strong molecular bands show rather weak dependence 

2 In fact these are spectral energy distributions (SEDs). 

3 For A > 2.7 /im model spectra shown here and in subse- 
quent plots were computed without molecular absorption; i.e, 
they are continuum fluxes, which provide information about 
the dependence of the continuum fluxes on input parameters. 
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Molecule D b c*l.E2 d*l.E6 e*l.E10 /*1.E15 



G 2 


6.116 


48.75 


.2192 


.4149 


.4121 


1.550 


CN 


7.700 


47.45 


.1332 


.1989 


.1778 


.6323 


CO 


11.105 


49.45 


.1651 


.3103 


.3100 


1.168 




Fig. 1. Model spectra of species that produced the 
strongest absoption features in the 1.0-2.5 ^m spectrum 
of Sakurai's Object spectra during 1997, computed for 
T e g /log g = 5500/0.0 model atmosphere with Asplund 
at al. (1999) abundances for October 1996. The model 
spectrum due to atomic species alone, labelled VALD (see 
text) is also shown. Spectra are artificially shifted on the 
y-axis. 

on log g, whereas the fluxes at 1.25-1.35, 1.60-1.75, 1.9-2.2 
microns show more noticeable changes. 

Previous abundance analyses of the spectra of 
Sakurai's Object and related R CrB stars indicate micro- 
turbulent velocities V t in the range 5-8 km/s (cf. Asplund 
et al. 2000). The value of Vt affects the spectral distri- 
bution, as is shown in Fig. ^. The effect of Vt on the IR 
spectra of Sakurai's Object is larger at the heads of molec- 
ular bands than elsewhere, because the heads are formed 
by closely packed molecular lines whose overall absorption 
is sensitive to V t . 

The main sources of line opacity in the model atmo- 
spheres approximating Sakurai's Object in 1997 are molec- 
ular (Pavlenko et al. 2000). Thus it is not surprising that 
the optical spectra which match Sakurai's Object respond 
weakly to changes in the hydrogen abundance. This is in 
contrast to the behavior of models corresponding to the 
star a year earlier (Asplund et al. 1997). Similarly, the 
model IR spectra of Sakurai's Object for T e g = 5000- 
6000 K depend weakly on log N(H) (Fig. |). The mag- 



Fig. 2. Dependence of the model IR spectrum on T e s . 
The model spectra use Asplund at al. (1999) abundances 
for October 1996. The observed spectrum of Sakurai's 
Object on July 13, 1997 is shifted on the y-axis. 
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Fig. 3. Dependence of the model spectrum on log g. 



nitude of the change in the spectrum when log N(H) is 
changed from -2.42 (the Asplund et al. 1999 value for 
October 1996) to -0.97 (i.e, a change of 1.5 dex) is com- 
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Fig. 4. Dependence of the model spectrum on V t . 



V4334 Sgr on April 21 
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Fig. 5. Dependence of the model spectrum of Sakurai's 
Object on log N(H) 

parable (in a qualitative sense) to lowering log g from 1 
to (Fig. |). 

3.3. Changes between 1997 April 21 and July 13 

Fits to the spectra of Sakurai's Object on April 21 and 
July 13 are shown in Figs. || and [7| The long wavelength 
portion of the H band is of special interest for the "car- 
bon problem," because the strongest absorption bands of 
the C2 molecule, the Ballick-Ramsey bands, occur just 
longward of 1.768 /im. In the comparatively hot atmo- 
sphere of Sakurai's Object log N(C) > log N(O) (Asplund 
et al. 1997, 1999) and the abundance of C2 depends mainly 
on the elemental abundance of carbon. Therefore, these 
bands may provide the most accurate determination of 
log N(C). The fits imply that the carbon abundance is in 



the range log N(C) = -1.3 ±0.2. The most likely value is 
0.3 dex higher than that found by Asplund et al. (1999). 
The accuracy of the determination of log N(C) is limited 
mainly by the quality of the molecular line list. 

The effective temperatures that best fit the 1.0-2.0 /im 
spectra in 1997 April and July are 5500 ± 200 K and 
5250 ± 200 K, respectively, indicating that the cooling 
evidenced by the dramatic spectral changes seen between 

1996 and 1997 (e.g., Geballe et al. 2002) continued in 1997. 
Our estimated uncertainties in the above temperatures are 
rather large, despite the comparatively good fits to the 
observed spectra, because of questions concerning abun- 
dances, non-sphericity effects, and dynamical phenomena, 
and because of contamination of the spectra by dust emis- 
sion (see below). 

3.4. Hot dust 

Emission by dust is evident in the 1997 spectra by the mis- 
match between the synthetic and observed spectra long- 
ward of 2.0 /im in Figs. ^ and ^. The difference between 
the observed and synthetic spectra is greater in the July 
spectrum, attesting to an increase in the amount of dust. 
The thermal emission from the dusty envelope overlaps 
the region of first overtone bands of 12 CO and 13 CO at 
A > 2.3 /im. Usually these bands are used for determina- 
tion of carbon abundances and isotopic ratios (cf. Lazarro 
et al. 1991). The reduced equivalent widths of the CO 
bands in July 1997 cannot be reasonably attributed to a 
large decrease in the oxygen abundance, because (1) this 
is unlikely to have occurred in three months and because 
the continuum shortward of the CO bands also shows an 
excess. We note that in fitting spectra, the most frequent 
situation is that the computed spectra have excess flux 
due to the deficit of known or hypothetical opacities. To 
fit the observed spectra, opacities in the model would need 
to be decreased at A > 2 /xm, an unrealistic possibility. 

4. Discussion 

Our analysis of the 1997 infrared spectra of Sakurai's 
Object strongly implies that dust was already present at 
that time. We note that Duerbeck (2002) did not find ev- 
idence for dust in the optical spectra from 1997. On the 
other hand, the fits by Pavlenko & Duerbeck (2001) to 
the observed SEDs at optical wavelengths indicate that 
Eb-v had increased by 0.6 (from 0.7 to 1.3) from April 

1997 to August 1998. However, the August 1998 data were 
best fit using T c g = 5250 ± 200 K, the same value as for 
July 1997 in this paper. Between June 1997 and August 

1998 there were some variations in the photospheric radi- 
aton, probably caused by mass losses events, evolution of 
the dusty envelope, and dynamical processes in the pho- 
tosphere -envelope system (see light curve of Sakurai's 
Object in Duerbeck 2002). Nevertheless, T c ff apparently 
remained nearly constant during this period. 

In 1997, the year of maximum optical brightness of 
Sakurai's Object, the luminosity then was still dominated 
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Fig. 6. Top: fits to observed spectrum of of Sakurai's 
Object on 1997 April 21. Bottom: details of the fits at 
1.6-2.0 /im; much of the structure at 1.82-1.95 /mi in the 
observed spectrum is due to incomplete removal of telluric 
absorption features. Synthetic spectra were computed for 
a microturbulcnt velocity of 6 km/s. 



by optical radiation. At that time "quasi-periodic fluctua- 
tions of increasing cycle length and amplitude were super- 
imposed on the general brightness evolution" (Duerbeck 
2002). In general, the effective temperature during such 
fluctuations does not need to follow changes of luminosity. 
In fact, it can be anti-correlated, since an increased radius 
can more than compensate for a lower T ff . On the other 
hand, a change of radius can change the thermodynamical 
properties in the radiating region (i.e. photosphere). That 
may explain the similarity of T c ff obtained in this paper 
for July 1997 and that found in August 1998 by Pavlenko 
& Duerbeck (2001). The decreased optical brightness of 
Sakurai's Object in 1998 was mainly caused by develop- 
ment of the dust envelope (Kimeswenger 1999). 



V4334 Sgr on July 13 
5250/0, log N(C) = -1.62 
log N(C) = -1.05 
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Fig. 7. Top: fits to observed spectrum of of Sakurai's 
Object on 1997 July 13. Bottom: details of the fit at 1.6- 
2.0 /an. Synthetic spectra were computed for a microtur- 
bulent velocity of 6 km/s. 



One question arises — were the optical and 2 /xm 
SED's being affected by the same dust in 1997-1998? The 
answer is probably yes. As mentioned earlier, the effective 
temperature remained constant during this time and thus 
cannot be the cause of the large change in E B -v • This 
suggests that the cause of the increase in E B _v was newly 
formed dust. The new dust would be expected to have 
been close to the star and thus quite hot. Indeed the full 1- 
5 /an spectrum from 1998 (e.g., Geballe et al. 2002) shows 
that the excess peaked close to 3 /im, indicating a mean 
dust temperature close to 1,000 K at that time. The dust 
must have been hotter (and closer) in 1997; this is sup- 
ported by the data from 1997 (Eyres et al. 1998; Geballe 
et al. 2002) which show that the continuum flux density 
decreased monotonically with wavelength in the observed 
wavelength range, 1-4 /im, at that time. Comparison of 
the 1997 and 1998 SEDs also show that much less dust 
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was present in 1997. Thus we confirm that the first ap- 
pearance of dust occured in 1997 and that the amount of 
dust increased through summer 1998, prior to its becom- 
ing totally dominant in the latter part of 1998 and since 
then. 
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